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ABSTRACT 

M87 is the first detected non-blazar extragalactic Tera-Electron-Volt (TeV) source with rapid variation and 
, very flat spectrum in the TeV band. To explain the two-peaks in the spectral energy distribution (SED) of the 

nucleus of M87 which is similar to those of blazars, the most commonly adopted models are the synchrotron 
' self-Compton scattering (SSC) models and the external inverse Compton scattering (EIC) models. Considering 

D , that there is no correlated variation in the soft band (from radio to X-ray) matching the TeV variation, and 

the TeV sources should not suffer from the 77 absorption due to the flat TeV spectrum, the EIC models are 
advantageous in modeling the TeV emission from M87. In this paper, we propose a self-consistent EIC model 
to explain the flat TeV spectrum of M87 within the framework of fully general relativity, where the background 
soft photons are from the advection-dominated accretion flow (ADAF) around the central black hole, and the 
high energy electrons are from the mini-jets which are powered by the magnetic reconnection in the main jet 
I't'I ', jGiannios et alJl20Toh . In our model, both the TeV flares observed in the years of 2005 and 2008 could be well 

l>H ■ explained: the 77 absorption for TeV photons is very low, even inside the region very close to the black hole 

' 20Rg ^ 50Rg; at the same region, the average EIC cooling time 10^ ^ 10^,?) is short, which is consistent with 

(~| , the observed time scale of TeV variation. Furthermore, we also discuss the possibility that the accompanying 

0^- X-ray flare in 2008 is due to the direct synchrotron radiation of the mini-jets. 

Q ' Subject headings: black hole physics — galaxies: individual (M87) 



1. INTRODUCTION 



Unlike blazars and the Galactic sources (e.g. SNRs), which composite the most population of TeV sources, M87 is the 
first discovered radio galaxy wit h TeV radiation. Recently three othe r radio galaxies Cen taurus A, 3C66B and ICS 10 were 
also identified as TeV sources dAharonian et al.ll2009t : lAhu et al.l 120091 : lAleksic et alj |2010). Comparing to blazars, M87 has 
much milder variations in both optical and X-ray band; most importantly, the prorninent kpc-scale j et of M87 has a very large 
T-j. viewing angle ~ 30° with respect to the line of sight (Bick nell & BegelmanllI996tlAharonian et al. | [2006). In 2005, rapid TeV 
, varia tion ( 1 - 2days) was disc overed with flat spectrum for the first time, but without correlated X-ray variation from the nuclear 
■ core dAharonian et al. Il2006h . While in 2008, the radio. X-ray, TeV joint observation discovered that TeV flares lasts 2 weeks 
• accompanied with an X-ray flare and a radio flare inside the unresolved nuclear core (30 x 6QRs) as well as with a radio blob 
moving out of the unresolved core (Acciari et al. 2009). Here used in the VLBA observations are based on the black hole 
(BH) mass of M87 Mbh = 6 x IO'^Mq (Gebhardt & Thomas 2009). Recently, a TeV flare in 2010 which is very similar to the 
previous TeV flares in 2005 and 2008 was reported. Further observations at X-rays and radio find the correlated X-ray flare, but 
no enhanced radio flu x from the core dAbramowski et al1l20I U [Harris et al.ll201 lb . In this paper, we adopt Mbh = 3 x IO^Mq 
dMacchetto et al .' 1997) to be consistent with the advection-dominated accretion flow (ADAF) models and the mini-jet models in 
the previous studies (Li et al. 2009; Giannios et al. 2010l). 
^ Although M87 has a large viewing angle, it is ver y near to the Earth with a distance of ^ 16.7Mpc. Owing to its proximity, it 
is proposed that M87 might be a misaligned blazar dTsvetanov et alJ[l998l) . Whereas blazars are believed to have jets beaming 
towards us, TeV sources in M87 could be shocks with relativistic bulk veloc ity inside the jets. Re cent observed minute-scale 
variation from galaxies like Mrk 421 (Fossati et al. 2008) and PKS 2155-304 dAharonian et al.ll200 7') indicates that TeV sources 
should be compact and very close to the BH. Comparing to the TeV flares of blazars, the slower variation (1-2 day) detected in 
M87 could be the results of too few observed TeV photons (each data point requires integration of photons for the whole night) 
or the much lower Doppler factor of the bulk velocity of the TeV sources due to the large inclination. Furthermore, unlike the 
steeper TeV spectrum of blazars, the much flatter TeV spectrum of M87 could be mainly due to the lack of 77 absorption. The 
large viewing angle may play an important role in avoiding 77 absorption within the jet; thanks to its proximity and the dimness 
of its host galaxy, the absorption to TeV photon from M8 7 caused by the galactic and th e intergalactic background soft photons, 
as well as the cosmic background photons is rather weak (INeronov & Aharonian 1120071) . 

To explain the spectral energy distribution (SED) of th e nucleus of M 87, the one zone syjichrotron self-Compton(SSC) models, 
which have been applied to the TeV flares in blazars (e.g. lTavecchio et al. 1998; Katarzyhski et al. 2001), face difficulties in fitting 
the two peaks in the SED, however, in the multi-zone SSC models, the TeV source is separated from the soft ones, therefore the 
whole SED c ould be easily fitted as long as there are enough SSC blobs with specific locations, electron distributions and bulk 
velocity (e.g. iTavecchio & GhiselUnill2008l: lLenainll2007t iGeorganopoulos et al.ll200 5V The multi-zone SSC models could also 
provide better explanation for the rapid variation, and the orphan TeV flares (the variations in the soft band can not connect with 
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Fig. 1 . — The contour of r^-^ = 0.1 or 1 in the q, /3 plane. Where r represents the distance to the BH in the Boyer-Lindquist frame. The viewing angle is taken 
to be 30° . The blank region in the lower panels is the region with r-,-, < 1 . 

the TeV variation). Furthermore, considering the very flat power-law SED of M87 in the TeV band with index p ~ -2.2 ~ —2.6, 
the 77 absorption provides crucial constrains on the SSC models. So far, most of the SSC models still have problems in explaining 
the very flat TeV spectrum due to th e certain 77 absorption, therefore the external in verse Compton (EIC) process could be more 
likely responsible for the TeV flare ("Begelm an etani2008t lOiannios et al.ll2009l 1201 0). In the EIC models, the TeV source is far 
away from the soft one, or has the relativistic bulk velocity with respect to the soft one, thus the 77 absorption can be reduced. 

In those models mentioned above, how the very high energy (VHE) particles are produced is still an open question. It is 
generally believed that VHE particles are accelerated by the shocks in jets. There are some other possibilities include the mini- 
jets powered via the magnetic reco nnection in the main jet fc iannios et al. 201 0i and referenc es therein); the magnetic centrifugal 
acceleration in the vicinity of BH (i Neronov & Aharonian ll2007i : iRieger & Agaronian ll2008h . and so on. In the mini-jets model 
(jOiannios et al. 2010), magnetic reconnection within the main jet produces two oppositely directed (in the rest frame of the 
main jet) mini-jets, in the laboratory frame, one of them always points within the angle of the main jet and is observable in 
blazars because their jets point at us. The other mini-jet (its counterpart) points outside the opening angle of the main jet and is 
potentially observable to off-axis observers in case of the misaligned jets, such as those of M87 and Centaurus A. In a word, the 
great advantage of the mini -jets model is that even at large inclination, mini -jets with a high bulk speed can still be detected, which 
is helpful for explanation of the fast TeV variation from M87. After all, in order to avoid the 77 absorption, the energy density 
of the soft synchrotron photons in the TeV source must be limited, therefore, the minimal Lorentz factor of the VHE particles is 
generally assumed to be 10^"^ and the strength of the magnetic field in the TeV source below several Gauss (G iannios et al. 2010) . 
Beside the direct inverse Compton process of the VHE electrons, there are also some alternative models to produce the SED of 
the TeV flares, such as the hadronic models including, among which, the interaction between the VHE protons and soft photons 
jReimer et alj|2004l) and the proton-proton collision process when a red giant was passing the base of the jet (Barkov et al. 2010). 

In this paper, we develop a fully general relativity EIC model for explaining the TeV emission from M87, in which the soft 
photons emanate from the ADAF around the BH. In §2, we investigate the safe zone of TeV photons (the 77 optical depth is 
below unity). The technical details of our fully general relativity EIC model can be found in §3. The numerical results and the 
discussions are given in §4 and §5, respectively. 

2. CONSTRAINTS ON THE LOCATION OF TEV SOURCES 

The rapid TeV variation of M87 (fvar < 1 - 2 days) found both in 2005, 2008 and 2010 indicates that the TeV photons should 
come from a compact source very close to the black hole (/?TeV ^ cfvarf^o, where (5d is the Doppler factor of the source). Especially 
the TeV flare in 2008, which is accompanied with a radio flare, lies inside the unresolved region and there is a radio blob moving 
outwards. Aga in, this indicates that the TeV source might be near to the BH and its distance might be less than ^ 100/?s 
(iBeiUcke et alJl2010h . 

As mentioned above, beside the constrains from the VLBA radio image and the variation time scale of the TeV emission, the 
background soft photons could also lay strong constrains on the location of the TeV sources, in consideration of the 77 absorption. 
Both the TeV flares in 2005, 2008 and 2010 had shown very flat power-law spectrum (0.1 ^ 10 TeV) strongly suggesting quite 
weak 77 absorption. Therefore, the EIC process see ms to be a more pl ausible mechanism to produce the flat TeV spectrum, while 
the SSC models suffer from certain 77 absorption jGiannios et al.ll2010i) . Even in the EIC models, if the source of soft photons 
is from a homogeneously isotropic blob whose size is about 50/Jg and infrared luminosity is about Lir ^ lO'^'erg/s, the location 
of the TeV sources could be still limited to be no deeper than ^ 5/?g inside that blob. 

The observed soft photons from M87 could be mainly either from the disk or from the outflow/jet. According to the VLBI 
observations of M87, its jet is well know n with a large incl i nation. However, the base of the jet seems to have a large opening angle 
(6'open ~ 30° inside the region 707?^, iJunor et alJl 19991: iLv et alJlIOOTt lAcciari et alll2009l) . therefore the jet might contribute 
to the observed soft emission. Due to the slight variation of the observed flux of soft photons (from radio to X-ray) from M87 
nucleus, here we suggest that the observed soft photons are mainly from the accretion disk around the central black hole as in Li 
et al. (2009). 
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By applying the ADAF model to the nuclear emission of M87, Li et al. (2009) obtained the accretion rate of ADAF in M87, 
generally consistent with the previous estimate of the Bondi accretion rate from the Chandra X-ray observation (Di Matteo et 
al. 2003). They then calculated optical depth of the radiation fields from the ADAF to TeV photons due to 77 absorption. The 
resultant optical depth suggests that the location of TeV sources sho uld be larger than ^ IQRg, in order to avoid the significant 
77 absorption to a 10 TeV photon (Wang et al. 2008; Li et al. 2009). Here we try to map a more detailed safe zone {t^j < 1) 
for TeV photons in the vicinity of BH. As in Li et al. (2009). For this purpose, we calculate the optical depth (r-^^) to 10 TeV 
photons emanating from vicinity of the central BH in the Kerr spacetime. The TeV photons move outwards along their geodesic 
trajectories and reach the observer's sky at points described by the impact parameters a and /3. Here a and (3 respectively represent 
the displacement of the image perpendicular to the projection of the rotation of the black hole on the sky and the displacement 
parallel to the projection of the axis (see e.g. Fig. 1 of Li et al. 2009). The resultant contour of t^^ as a function of the location 
of the TeV sources for r-y^ = 1 and 0.1 in the a/S plane can be obtained (see Fig.fTli. 

As shown in Fig.[Tl it is obvious that the larger the spin a, the deeper the safe zone. For instance, for a = 0.998, the safe zone 
can reach even to IQRg along some trajectories of lOTeV photons, the reason is that for larger BH spin, the infrared-UV radiation 
is concentrated in the inner disk. Therefore, the colliding angle between the soft photons and the TeV photons in the mini-jets is 
smaller, significantly reducing the 77 absorption. Besides, in the panels of a = 0,0.8, the asymmetricity is caused by the rotation 
of the disk and the viewing angle(30°) of the observer While in the panels of a = 0.998, since the disk become a more compact 
source of soft photons, the asymmetricity caused by the disk become weaker 

3. FULLY GENERAL RELATIVITY EXTERNAL-INVERSE COMPTON MODEL 
3.1. The TeV source 

In our EIC model, following iGiannios et al.l (120101) . the mini -jets powered by magnetic energy in the jet correspond to be the 
TeV sources. For a clarity, we present the details of the model in what follows. 
The strength of the magnetic field in the main jet is estimated as. 



R / r 



4tt [A-nrjcT^j I \ l + (J 



(1) 



where Bj is the magnetic strength, Lj iso (~ 10'*^ erg/s) is the observed isotropic power of the jet of M87 (iBicknell & Begeknan I 
1 1996t [Reynolds et al. lll996l:IOwen et al.ll2000l: IStawarz et al.ll2006t iBromberg & Levinsonll2009^ ■ Aj (Trrj^) is the cross sectional 
area of the jet, a is the ratio of the magnetization energy to the kinetic energy of the main jet, Fy is the Lorentz factor of the main 
jet, and c is the light speed. 

Although the inclination of the jet is about 30", the mini-jets (blobs of plasma with characteristic Lorentz factor Tco at angle 
^mini with rcspect to the jet in the jet rest frame) are capable of beaming towards us with a high bulk velocity Fem = FjFco(l + 
^j/3cocos0mini) in the laboratory frame. The bulk Lorentz factor of the mini-jets in the jet rest frame is about T^o ^ cr'/^, which 
corresponds to the Alfven speed of the plasma in the jet. The characteristic thermal Lorenz factor of the electrons in the mini-jets 
rest frame is about ^ fa^l^m^/m^. The energy distribution of the high energy electrons is assumed to be 

K{i) = Ne{n)/5l=N^,iP, (2) 

where 10"^ "^l' = 1 1 So < 00, p < -3.2. In our canonical model, Fy is taken to be 5, and Fr^ = 12 (Sn = 23, if the mini-jets are 
beaming toward us). These parameters are consistent with the original mini-jets model of IGiannios et al.l ( 120101) . 

3.2. The Local Soft Radiation Field 

In our EIC model, we assume that the soft photons are from the accretion disk and the TeV source (mini-jets) is located along 
the major axis of jet with a height tfrev ab ove the central bla ck hole. To investigate the EIC process, we apply the same ray-tracing 
technique as discussed in Li et al.l (l2009h . In iLi et al] (12009 ). the authors discussed the optical depth of the TeV photons due to 
colliding with the soft photons from the accretion disk. In this work, we investigate the Compton scattering of the soft photons 
from the disk by the VH E electrons in the mini-jets (as shown in Fig.|2]i, therefore, we just simply replace the TeV photons in the 
model ofl Li et all (120091) with the VHE electrons. 

To obtain the flux density of the soft photons from the disk, the global dynamical structure of the ADAF, such as the four 
velocity of the fluid in the disk, the temperature of ions and electrons, should be determined first; then the local emergent spectra 
Iiy,Md) at the radius in the rest frame of the fluid can be calculated. Using the ray-tracing technique, the obse rved SED can b e 
obtained to fit the multi-wavelength observations of M87. The fitting parameters of the ADAF are obtained by lLi et aP (l2009l) . 
In what follows, we summarize the procedure to obtain the radiation energy density of the soft photons from the direction (0s, 0s) 
at the colliding location {rc,dc,(pc)- 

In the locally non-rotating frame (LNRF), at the interacting place (rc,9c,4>c) , the two motion constants of the soft photons (As, 
Qs) are related with their traveling direction (9^, (/"s) as follows: 

As = — Qs = ^'-(flCOS0c)' + (AsCOt0e)'. (3) 

where, 

sin 6, sin (h ^ sin 9 c A ^ sin9,cos6A^^-^(l-(jjX,) 
£/ = — — <i§ = — (4) 
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Fig. 2. — A schematic picture for the disk dominating external Compton-scattering model. A mini-jet is beaming toward the infinity observer, and the soft 
photons are from the disk. 



where A, S, A, and lo are the metric functions defined in lBardeen et alJ (Il972h . 

After knowing the consta nts of motion As and Qs, the soft photo ns can be traced back to the disk at certain radius by solving 
the geodesic equations (e.g. iBardeen et al.lll972t lYuan et ani2009l) 



sr = ± 



dr 



= ± 



(5) 



where =5^ is the affine parameter and the ± signs represent the increment (+) or decrement (-) of r and 6 coordinates along the 
trajectory, respectively. 

Along these trajectories, the redshift factor for a soft photon travels from disk to the IC interaction location can be obtained 
by 



e(';,(LNRF)P;, 



e{;)(LRF)P^ 



70^ 



1/2 

17^ 



(6) 



where and are the frequency of the soft photon at the colliding place and the disk, and 7,. = (1-/3^) '/^ and 7^ = (1 -Z^^) '/^ 
are the Lorentz factors of the radial and azimuthal velocity of the fluid in the accretion disk, respectively. 

According to the Liouvell theorem, the final radiation energy density of the soft photons at the IC location (rc,6'c,0c) can be 
written as 

/i.j(/Wd,'-d)^s 



(7) 



where h is the Planck's constant. 



3.3. The VHE electrons 



After determining the radiation energy density of the soft photons Us at the colliding points, the direction of motion of the 
VHE electrons which produce the observed TeV photons is needed to calculate the TeV spectra. Due to the effects of special 
relativity, it is a reasonable assumption that both the direction of the relativistic electrons and that of the TeV photons are the 
same. Denoting the two constants of the motion of the TeV photons as Aic and Qic, the beaming direction of the TeV photons is 
as follows: 
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where the trajectory of the TeV photons determined by Aic and 2ic can be found by tracing the observed TeV photons from the 
infinity (r = cx),f?obs = 30°) to the interacting location (r = r^^O = 9^). As assumed above, the TeV source is located inside the jet, 
therefore, dc = Q and subsequently Aic = 0. 

After determining the trajectory (Aic,2ic) of the TeV photons, the corresponding redshift factor gic for the 7-ray photons 
traveling from the IC interaction location to the infinity is given by 
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3.4. Spectrum of Compton-scattered external radiation fields 
Given the radiation ener gy of the soft pho tons and the energy dist ribution of the high energy electrons, the Compton spectral 



Uiven tne raaiation ener gy 01 tne sort pno tons ana tne energy aist ri button or tne nig 
luminosity is given by (see lFinke et ani2008l and iDermer et al.lE009l for more details): 



fhcai(<^ic) = 
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where 65 is the energy of the soft photons, ejc is the energy of the 7-ray photons created via the IC process, /^s = cos(6's), and 
flic = cos(6'ic). With the approximation that scattered photons travel in the same direction as the VHE electrons, the Compton 
cross section can be drawn as (Dermer et al. 1993 ; Dermer et al. 2006) 
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e = 7es(l- 
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and V'lc is the interaction angle between the VHE electron and the soft photon during the IC process. 



cos V'lc = /^ic/^s + V 1 ~ M/c y'^~K cos(0/c - 0.,) 



(13) 
(14) 

(15) 
(16) 

(17) 



The optical depth T(eic,//Tev,fl) of these 7-ray photons (eic) can be integrated along their trajectories as in Li et al. (2009), 

7"77(Aic,2ic,eic) = JJJ (l-cos^^^)ajj(es,eic,i^j'f)^8ldflsdesdl, (18) 

where dl = e'^Ed.T is the proper length differential with d,3'^ defined to be differential of the affine parameter ,3^ along the 
trajectory of the TeV photons, and tp^^ is the interacting angle between the TeV photon (ejc) and the soft photon (es). 
Consequently, the final observed SED in the infinity would be, 

/IS* (^ic^) = [{f-^'fRU^iO. (19) 

4. NUMERICAL RESULTS 

4. 1 . //xev , fl dependence 

In this paper, we choose the BH spin as a = 0, 0.8, 0.998 to be consistent with the ADAF model used in fitting the SED (from 
radio to X-ray) of M87 in Li et al. (2009). The location of the TeV source is set to be very close to the BH as //xev = 5, 10,20,50, 
and 100/?^, respectively. The distribution of the VHE electrons is described by A^o = 0.4 x 10^", 7 ^;^ = 10"^, and p = - 3.2, where 
the maximal index p is chosen to fit the very flat TeV flare spectrum observed in 2005 and 2008 tGiannios et al.ll20Toh . 

In Fig. [3] we show the dependence of the SED of the EIC scattering on the height of the TeV source and BH spins. We can find 
that the efficiency of the Compton-scattering increases dramatically if the TeV source is closer to the disk. Furthermore, if the 
location of the TeV source //xev is below 2QRg, the Compton-scattered luminosity of the soft radiation fields from the the disk 
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Fig. 3. — SED of the external Compton-scattering model with the different spins of the black hole (a = 0,0.8.0.998) and the location of the mini-jets is above 
the black hole {Left panel: //xeV = 5, 50/?g, right panel: //tcV = 10, 20, lOORg). The energy distribution of the VHE electrons are determined by thi'ee pai'ameters: 
Nq = 0.4 X lO'", 7min = 10** and p = —3.2. The Lorentz factor of the mini-jets in the laboratory frame is taken to be Fcm = 12. MAGIC observation in 2008 is 
shown in the symbol of butterfly (Albert et al. 2008). 

around the black hole with spin a = 0.998 is higher than that from the black hole with spin a = 0; meanwhile, if the location of 
the TeV source is higher than 2QR^, for a = 0, the Compton-scattered luminosity is higher than that for a = 0.998. This is because 
for a = 0.998, the soft radiation field from the disk is more compact, i.e., the soft radiation field is strong in the vicinity of the 
BH, however, at the distance far away from the disk, the region of the soft radiation field looks like a point source. Therefore, the 
colliding angle between the soft photons and the VHE electrons in the mini-jets is smaller, which reduces the rate of the inverse 
Compton scattering significantly. 

It is also clearly shown in Fig.[3]that each spectrum has an exponential cut off beyond Ecut ^ lO^^^mgC^ which is caused by the 
77 absorption. If the location of the TeV source is nearer to the BH (//xev < 20/?g), the spectrum will also suffer from the severe 
77 absorption in the observable band 0.1 - lOTeV. In each of those reduced spectrums of //jev = 5, 10, there is always a warp up 
trend at Eic ~ lO^WeC^ which is caused by the relatively low optical depth t^^ at that band, since the spectrum of the disk has a 
relatively low luminosity at ^ lO'^mgC^. 

4.2. TeV flares and The X-ray companion 

In 2008, observations had shown that the TeV, X-ray, and radio flares are well correlated (l Acciari et al.ll200 9*). though no X-ray 
correlation was found in 2005. We propose that the accompanying X-ray flare in 2008 results from the synchrotron radiation of 
the VHE electrons in the mini-jets, in other words, the X-ray flare could be a byproduct of the TeV flare, which just happens 
to be above the background X-ray flux. While in 2005, it is more likely the X-ray flaring flux is below the more stable X-ray 
background which is attributed to the disk and/or the corona. It is notified that the TeV flare in 2008 is stronger than that in 2005, 
and the corresponding X-ray flare is rather mild: during the flaring, the X-ray flux only increased about two times above the 
average. 

In Fig.m both the SED of the direct synchrotron radiation and the Compton-scattered radiation are shown. We choose //xev = 
20,50 to be as close to the BH as possible while avoid severe 77 absorption. Owing to the weak 77 absorption, all the EIC 
spectra have inherited the original power-law index of the VHE electrons, which is p = -3.2,-3.5. As expected, the corresponding 
synchrotron spectrum fitting the X-ray flares in 2008 is comparable with or exceeds the average X-ray flux, while the X-ray flux 
in 2005 is buried beneath the background flux. 

The magnetic field in the mini-jets is taken to be B = 1.6,4,2.6,8G, respectively, which is used to calculate the direct syn- 
chrotron radiation. The strength of the magnetic field used in this work is consistent with that in the original mini-jet model 
Bern (iGiannios etal.l i2010. and references therein), that is, < 0.8(Lj,iso/10'*^ erg/s)'/^(100j;p/ rjet)(5/r;)^ G. We propose that the 
base of jet is cone shaped as rjet ^ O.SHsindopen, where the full opening angle ^open ^ 30° (lAcciari et al.ll2009l) . which leads to 
B{50Rg) < 6.4G, B(20Rg) < 16G. The more detailed parameters are shown in Table 1. 

There are four main parameters NQ,p,Hi^\j ,a in our models. Where the spin a is fixed because it only slightly influence the 
SED for //xeV = 20,50 as shown in Fig.|3j the electron power-law index p is fixed to follow the observed TeV spectrum index 
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Fig. 4. — Synchrotron (dashed Hnes and dash-dotted hnes) and the external Compton-scattered (solid lines and dotted lines) emission from a Mini-jet which are 
modeled to fit the observed emission from M87, the detailed model parameters are shown in Table one. The observed emission include: (1). MAGIC observation 
in 2008 which is shown in the symbol of butterfly (Albert et al. 2008); (2). HESS observation in 2005 which is indicated by the black square with long error-bars 
; (3). HESS observation in 2004 which is labeled as the black triangle with long en'or-bar (Aharonian et al. 2006); (4). FERMI observation in 2009 which is 
labeled as the black solid circle with error-bars (Abdo et al.2009); (5). The approximate data from Chandra is shown in the hollow triangles. With 5% uncertainty, 
the one on the top repres ent the flux of nucleus in 2008 Feb 16 (data are provided by Dr. D. E. Hanis) the other three rep r esent the flux in 2005 Apr 22, Apr 
28, May 04 se parately I H ams et alj|2006l , l2009[) ; (6). Typical nucleus emission is shown in the asterisk. fHarms et al. lll994l ; IMaoz et al. 1120051 ; iRevnolds et al~l 
II 99ft [Sparks et al. 1996). The synchrotron-self Compton emission is too weak to be shown in this figure. 



TABLE 1 

Model parameters used in fitting the observations shown in Fig. 4. 



EIC hnes in Fig. 4 


Hj^yiRg) 7n,i„ A'o(counts) 


P 


B(G) 




Fitting 2008MAGIC'' with two lines on the top 






dash-3doted line 
dotted line 


50 10" 1.0 xlO^' 
20 lO'' 1.0 xlO'" 


-3.2 
-3.2 


2.6 
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Fitting 2005HESS'' with two lines in the middle 






dashed line 
dash-dotted line 


50 10'' 0.6 xlO^' 
20 10-* 0.4 X lO'" 


-3.2 
-3.2 


1.6 
4 




Fitting 2004HESS''+2009FERMI'' with two lines at the bottom 




solid line 
sohd line 


50 10"* 1.6x10" 
20 10" 0.5 X 10^' 


-3.5 
-3.5 





Note. — The bulk speed of the mini-jets in local lab frame is taken as T^i^i = 12 and the viewing angle is set as Sobs = 30°. 
The parameter spin of the BH is chosen as a = 0.8 , because the changing of the spin barely influence the shape of the spectrum, 
it just slightly boost the entire SED up or down when applying //tcV = 20, 50Sg. 

" MAGIC observation of M87 in 2008, as a flaring state 
HESS observation of M87 in 2005, as a flaring state 
HESS observation of M87 in 2004, as a quiet state 
FERMI observation of M87 in 2009, as a quiet state 

and the 77 absorption is too low to soften the spectrum; the other two parameters A^o and //xev are degenerated, making it hard 
to draw any soHd constrains by fitting the SED. To be more specific, the final TeV flux should be roughly Fic ^ No/Hj^y. As 
shown in Table 1, to fit the same observation data, when the i/xev changes from 207?g to 50/?g, we need about 10 times more VHE 
electrons. Therefore, considering that both the TeV flares in 2005 and 2008 have the same a and p{-2.2 ^ -2.6), we suggest that 
the TeV source in 2008 may be more powerful or/and be closer to the BH. 

Being closer to the BH, the TeV source will suffer more 77 absorption, which would soften the spectrum. But no obvious 
changes of the spectrum index can be drawn from the observations in 2005 and 2008 due to the long error-bars. It is noticed that, 
as shown in Fig. 3, the 77 absorption barely influences the 7-ray spectrum below 0.1 TeV, therefore the Fermi telescope should 
be able to provide better constrains on the power-law index of the VHE electrons. 

As shown in Fig. g] we also fit the quiet state (2004HESSH-2009FERMI) spectrum of M87 with the mini-jets model. In 
other words, the quiet TeV flux could be attributed from some mini-jets which are weaker, misaligned, or heavily absorbed by 
synchrotron or/and 77 pair production (lower Ne,p,5£,)- Where the misaligned Mini -jets could be considered as weaker mini- 
jets consisted of much fewer VHE electrons beaming towards us with lower bulk velocity, see details in (.Giannios et al...2010.) . 
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Therefor we could not exclude the possibility that weaker flares are due to misaligned mini-jets, which are supposed to have 
lower Sd. It should be noticed that lower Su could not explain the steeper spectrum. 

As mentioned above, the soft photons field are demanded to satisfy both the low 77 absorption rate and the high IC cooling 
efficiency. In models as shown in Fig. HI apparently, the synchrotron cooling effect is not dominating since the EIC flux is much 
higher. The IC cooling time is short enough so that there is no influence on the observed rapid variation (1 ^ 2days): 



where Lcev-Tev is the isotropic total EIC flux and is the total energy of VHE electrons in the lab frame. With the parameters 
listed in Table 1, as expected, one can find fcooi ^ 10^~^s. During this cooling time, the mini-jets can only travel less than ^ Q.5Rg, 
which is consistent with the above assumption that the location of the TeV source keeps unchanged during a flare event. 



In this paper, by taking into account the fully general relativity effects, we propose a disk-dominating external Compton- 
scattering model for explaining the flat TeV radiation from M87. The external Compton-scattering model suffers much less self- 
77 absorption, comparing to the one-zone self-synchrotron Compton model, in which the soft photons and the VHE electrons 
are from the same blobs. The advantage of the EIC model in which the soft radiation is from the accretion disk is that the 
soft radiation source could be more compact and further away from the TeV source than ones in the SSC models, in which 
the 77 absorption can be reduced significantly. Besides, our EIC model is also supported by the observed mild IR-UV-X flux 
variations of the nucleus and the large viewing angle of the jet of M87 (the nuclei flux with viewing angle 30° is unlikely to be 
contaminated by the jet, neither to be blocked by the disk itself). It should be noticed that the VHE electrons and TeV photons 
may be influenced by some outflows with the anisotropic radiation at the soft band which are not beaming towards us, such as 
the jet itself. Unfortunately, the base of the jet of M87 ap pears to be quite chaotic, thus it is hard to draw any conclusions about 
the detailed structure of jet/outflows jPerlman et al.ll201 lb . In our model, the TeV sources are the mini -jets which are easily able 
to beam towards us with high Doppler factor and give rise to the observed VHE radiation. 

How to distinguish between SSC and EIC models by the future observations? The main difference between these two models 
is that the IR-UV variation and TeV flares are not necessarily correlated in the EIC models, but they are in the SSC models, 
especially the multi-zone SSC model. Therefore, we can distinguish between the SSC and EIC models by the simultaneous 
observations of the IR-UV emission during the TeV flare. Unfortunately, during the 2005 TeV flares, there is no accompanied 
HST observation of M87. With the better sensitivity of CTA (it is about >10 times better than those of the present Cerenkov 
telescopes), we will be able to detect the fainter flares and obtain the minimal variation time scale of the TeV flares, which 
increases the opportunity of finding the correlated variation between IR-UV and TeV emission. 

In our models, with little 77 absorption, all the EIC spectra at the observational band have inherited the original power-law 
index of the VHE electrons. After fitting those spectra, we calculate the cooling time of the VHE electrons in the mini-jets (see 
section 4) to make sure that they would not travel longer than cty^So- It turns out when TeV source is located at //jev < 50/?g, 
the IC cooling time is very much shorter than the upper-limits of the observed TeV variation of 1-2 days. Therefore, if possibly, 
the TeV source is located further away from the disk (//jev > 5QRg), the cooling efficiency of the VHE electrons in the mini-jets 
is much lower, the mini-jets could move out of the main jet. As a result, the variation of the soft radiation background at the 
different locations above the disk should be considered. 

We also discuss the probability that the direct synchrotron radiation from the mini-jets may cause the correlated X-ray flare in 
2008. As expected, according to our model, if the magnetic field is about several Gauss, the direct synchrotron flux from these 
mini-jets which lies at the X-ray band can explain the X-ray flare very well but has no influence on the observed radiation at the 
GeV-TeV band. However, the magnetic field B can not be constrained by any direct observation, as discussed above, and we only 
have a rough estimate of the strength of the magnetic field supposed to be consistent with the mini-jets model. Considering that 
the synchrotron flux is roughly about F^yn ^ NqB^, it is obvious that the A^o and B are degenerated. If the TeV source is confined 
inside the region 20Rg < i/xev < 50/?g, the corresponding strength of the magnetic field is required to be about 2G < B < 8G, in 
order to produce the X-ray flare in 2008. 

Dealing with the two weeks data of the 2008 flare from MAGIC, Albert et al. (2008) had separated the high state data ( ^ 
1 day time scale) from the low state data and they found that during the high states, the spectra are harder: the spectrum index 
of the high state is about 2.2, while that of the low state is about 2.6. Because the IC cooling time and the synchrotron cooling 
time of VHE particles is short, the observed fastest variation should be dominated by the time of the acceleration, or even the 
shifting of the beaming angle of the mini-jets. If there is no continuous acceleration of the relativistic electrons, the spectra 
of the TeV flares will become steeper In the year of 2008, VERITAS captured the TeV flare (Feb 9-13, 2008) which is close 
to the date of the X-ray flare (Feb 16, 2008). VERITAS obser vation shows that the spectrum index of the TeV flare is about 
2.4, which is softer than the earlier flares captured by MAGIC ( lAcciari et al.l l2010). We suggest that the TeV flare captured by 
VERITAS is happening during the rapid synchrotron cooling of the relativistic electrons in the mini-jet, and the X-ray flare could 
be a byproduct of that TeV flare, as the X-ray flares in 2008. Our suggestion could naturally explain the X-ray excess and the 
softer spectrum of TeV emission. There are two ways to check above suggestion: first, the flux of the X-ray via synchrotron 
radiation should be anti-correlated with the TeV flux via EIC, which could be tested by the future observations with the higher 
time resolution (less than 1 day); second, and the TeV flare s hould share the sim ilar spectrum index with that of the X-ray flare. 
Unfortunately, due to the pile up effect of X-ray observation ( [Harris et al.ll2009l) , we can not obtain the correct spectrum index of 
the X-ray flare. 




(20) 



5. CONCLUSIONS AND DISCUSSIONS 
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